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Abstract 

The B c — > J/tpTT, r] c ir decays are studied with the perturbative QCD approach. It is found that 
the form factors and F^ c ~*^ c for the B c — > J ftp, 7] c transitions and the branching ratios are 

sensitive to the parameters uj, v, fj/^p and f Vc , where u and v are the parameters of the charmonium 
wave functions for Coulomb potential and harmonic oscillator potential, respectively, fj/^ and / r?c 
are the decay constants of the J/ip and r) c mesons, respectively. The large branching ratios and the 
clear signals of the final states make the B c — ► J/i/jtv, rj c ir decays to be the prospective channels 
for measurements at the hadron colliders. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) is schedule to run in this year. At the era of the LHC, 
there is still a room for B physics. The study of the decays of B mesons is important and 
interesting for the determination of the flavor parameters of the Standard Model (SM), the 
exploration of CP violation, the search of new physics beyond SM, and so on. The decays of 
B u d mesons have been investigated widely by the detectors at the e + e~ colliders, such as the 
CLEO, Babar, Belle. The B c meson could be produced abundantly and studied detailedly 
at the hadron colliders, such as the Tevatron and LHC. The study of the B c mesons will 
highlight the advantages of B physics. 

Compared with the B u & mesons, the B c mesons have some special properties: (1) The 
B c mesons are the "double heavy-flavored" binding systems. We can study the two heavy 
flavors of both b and c quarks simultaneously with the B c mesons. (2) The B c mesons 
have much rich decay modes, because they have sufficiently large mass and that either b 
or c quarks can decay individually. The potential decays of the B c mesons permit us to 
over-constrain quantities determined by the B U ( i meson decays. 

It is estimated that one could expect around 5 x 10 10 B c events per year at LHC Jl|]. The 
nonleptonic decays of the B c mesons have been studied in previous literature l[ Q]. The 
theoretical status of the B c meson was reviewed in l|. In this paper, we will concentrate 
on the B c — > J/ipir, t] c ir decays using the perturbative QCD approach. There are several 
reasons : 

(i) From the experimental point of view, the decay modes containing the signal of J/if> 
meson are among the most easily reconstructible B c decay modes, due to the narrow- 
peak of Jfij) and the high purity J/ip — > £ + £~ . For example, the B c mesons are firstly 



discovered via B c -> Jj^lv by the CDF Collaboration in 1998 [3]. Recently the CDF 
and DO Collaborations announced their accurate measurements on the B c mesons via 
B c — > J/ipir mode 0, 0]. Especially, compared with the semi-leptonic decays where 
the neutrino momentum is not detected directly, all final-state particles are detectable 
for the B c — > J/ipir, i] c tt decays. It is estimated that the ATLAS detector would be 
able to record about 5600 events of B c — > J/ipn per year l|. So B c -> J/ipir, r] c ir 
decays may be two of the most prospective channels for measurements. 
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(ii) From the phenomenological point of view: In recent years, several attractive methods 
lave been proposed to study the nonleptonic B decays, such as the QCD factorization 
(|, perturbative QCD method (pQCD) 0, [si, ^], soft and collinear effective theory 
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111 ], and so on. The study of B c decays provides opportunities to test the and 
collinear factorizations, to check the various treatments for the entanglement of differ- 
ent energy modes, to deepen our understanding on perturbative and nonperturbative 
contributions. These methods developed recently are widely applied to the nonleptonic 
two-body B u ^ s decays in literature, but with very few application of these methods 
on the B c meson decays. The appealing feature of the pQCD factorization 
is that form factors can be computed in terms of wave functions (nonperturbative 
contributions) and hard kernels (perturbative contributions arising from hard gluon 
exchange) assuming that additional soft contributions are suppressed by the Sudakov 
factor in the heavy quark limit. Although there is still some controversy about the 



pQCD method, for example, the problem of gauge invariant [12], the pQCD method 
has been extensively used in the past to study nonleptonic B decays with fairly good 
phenomenological results [3]. In this paper, we will take the B c — > J/ipir, i] c 7r decays 
as examples to discuss the B c decays in the perturbative QCD method. 

(iii) From the theoretical point of view: The B c — > J/ip7i, i] c tt decays are similar to the 
B q -> £>Wtt (where q = u, d, s) decays with the "spectator quark" ansatz. The 



<i 

B q — > D q *^n decays have been studied with the pQCD method 15]. Compared with 



the B q — > D^tt decays, the B c — ► J/4>tt, t] c tc decays are easy to deal with because 
that the B c meson and the J/ip (or i] c ) meson are heavy quarkonia and could be 
described approximatively by nonrelativistic dynamics. Given mg c ~ m& + m c , the 
wave function of the B c mesons would be close to 5(x — m c /mB c ) in the nonrelativistic 
limit (where ms c , mj, and m c are the masses of the B c mesons, b quark, and c quark, 
respectively; x denotes the momentum fraction of the c quark in the B c meson). The 
wave functions for pion are well-defined. The only parameter is the wave function of 
the J j ip (or r] c ) meson. So the B c — ► J/ ipn, r) c ir decays provide good platform to test 
quark potential models derived from QCD. 

This paper is organized as follows: In Section [Til we discuss the theoretical framework and 
compute the decay amplitudes for B c — > J/ipn, t] c tt with the perturbative QCD approach. 
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The section IIHI is devoted to the numerical results. Finally, we summarize in Section IIVI 



II. THEORETICAL FRAMEWORK AND THE DECAY AMPLITUDES 



A. The effective Hamiltonian 



Using the operator product expansion and renormalization group (RG) equation, the low 
energy effective Hamiltonian for B c — > X cg 7r decay can be written as (where X c5 = J/ip, i] c ): 

n eff = ^V A V^{C x {AQi + CMQ 2 } + H.c, (1) 

where V c bV* d is the CKM factor accounting for the strengths of the concerned nonleptonic 
decay processes. The parameters Cj(/i) are Wilson coefficients which have been evaluated 
to the next-to-leading order with the perturbation theory. The expressions of the local 
operators are 

Ql = [CaT/X 1 - 75)UK37 M (1 - 7 B )U/9], Q2 = [Calni^ ~ 75)MK?7 M (1 - 75^], (2) 

where a, j3 are color indices. The essential problem obstructing the calculation of decay 
amplitude is how to evaluate the hadronic matrix elements of the local operators. 



B. Hadronic matrix elements 



The calculation of the hadronic matrix elements is difficult due to the nonperturbative 
effects arising from the strong interactions. Phenomenologically, the simplest approach 



to hadronic matrix elements is the Bauer-Stech-Wirbel (BSW) model 14J based on color 
transparency and naive factorization hypothesis, where the hadronic matrix elements are 
parameterized into the product of the decay constants and the transition form factors. 
One defect of the rough BSW method is that the hadronic matrix elements cannot cancel 
the renormalization scheme- and scale- dependence of the Wilson coefficients. To remedy 
this problem, the "nonfactorizable" contributions must be taken into account. Using the 
Brodsky-Lepage approach [16], the hadronic matrix elements can be written as the convo- 
lution of a hard-scattering amplitude, including some perturbative QCD contributions, and 
meson wave functions. 
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Recently, a modified perturbative QCD formalism has been proposed under the kx fac- 
torization framework [t^S, [9]. The Sudakov effects are introduced to modify the endpoint 
behavior. The decay amplitudes are factorized into three convolution factors: the "harder" 
functions, the heavy quark decay subamplitudes, and the nonperturbative meson wave func- 
tions, which are characterized by the W ± boson mass myy, the typical scale t of the decay 
processes, and the hadronic scale Aqcd, respectively. Using the resummation technique and 
the RG treatment, the final decay amplitudes can be expressed as 

A(B c ^X c5 tc) cc C(t)®H(t)®$ B M, h)^x cS (x 2 , b 2 )^(x 3 , b 3 ), (3) 

where the Wilson coefficient C(t) is calculated in perturbative theory at scale of mw and 
evolved down to the typical scale t using the RG equations, ® denotes the convolution over 
parton kinematic variables, H(t) is the hard-scattering subamplitude, the wave functions 
b) absorb nonperturbative long-distance dynamics, x is the longitudinal momentum 
fraction of the valence quark of the meson, b is the conjugate variable of the transverse 
momentum of the valence quark of the meson. According the arguments in 

Q, fl 3, the 

amplitude of Eq. ([3]) is free from the renormalization scale dependence. 



C. Kinematic variables 



For convenience, the kinematics variables are described in the terms of the light cone 
coordinate. The momenta of the valence quarks and hadrons in the rest frame of the B c 
meson are defined by 

pi = =Jf (1,1,00, h = x lPl , n 2 = (1,0,0), 

P2 = (1, 4 eg , 00, k 2 = x 2 p 2 + (0, 0, k 2± ), ey = ^(1, -rj, 0), 

p 3 = (0, 1 - r 2 XcS , 00, h = x m + (0, 0, k 3± ), n 3 = (0, 1, 0), 

where the notation of momenta of pi and k{ are displayed in FIGiB The null vectors n 2 and 
n 3 are the plus and minus directions, respectively. The mass of the 7r meson is neglected. 
The momentum of the 7r meson is chosen to be parallel to the null vector n 3 . The mass 
ratios are r Xc _ = m Xc Jm Bc , r b = m b /m Bc , r c = m c /m Bc . 
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D. Bilinear operator matrix elements for mesons 



In terms of the notation in [17(, the nonlocal bilinear-quark operator matrix elements 
associated with the B meson, , meson, the longitndinaLy polarized J H meson, r k meson 
are decomposed into [17|, [18[ 



\c a (z)b p (0)\B;( Pl )) 



(n-(p 3 )\d a (0)u p (z)\0) 



Pot 



—I 



y/2N c J 

(J/^(p 2 ,e|,)|c a (0) C/3 (z)|0) 
(Vc(P2)\c a (0)cp(z)\0) 



dx 3 e + "3 p 3 * [ 75 [ fa<f>* (x 3 ) (x 3 ) - fa (ft 3 /i 2 - n r n 2 ) 0*(ar 3 )l } , 

i ^ L \ ) pa 



V2K 



Pa' 



(4) 

(5) 
(6) 

^2W C J ~ l^l^^v-/ '-vc^vwjj^. ( 7 ) 

where the wave functions 0", 0^, are twist-2, 0^, 0^, 0^,, 0* c are twist-3, /i^ = m^/(m u + 
m^) . Their expressions are collected in APPENDIX |A] and [B] 

For the wave function 0s c , we will take the nonrelativistic approximation as stated in the 
introduction, i.e. 

A IB, 



^ 2V2K 6{X rc) ' 
where N c is the color number, fs c is the decay constant of the B c meson. 



(8) 



E. B c — » X C c form factors 



The B r — > X r? form factors are defined as 



{r] c {p 2 )\^b\B c { Pl )) 
i i 



14 



19|: 



2 2 

mi, — 



(J/^(p 2 ,e)|c 7 ^75&|5 c (Pi)) 



(9) 



(10) 



where q = p± — p 2 , e* denotes the polarization vector of the meson. Fq[^ Vc and A B( {^ 3 
are the transition form factors. In addition, at large recoil limit, q 2 = 0, we have 

F(f c ^ c (0) = Ff^ Vc {0), A^ J/ ^{0) = Af-^(0), (11) 
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A^ J, \q 2 ) = mB ° +m W A f^Wtf) _ m B-™>J/j, A B c ^ {q2) (12) 
In the perturbative QCD approach, these form factors can be generally written as 

A B^J/iP (()r F B c -> Vc) ^ QB^bjQHWQQx^xM- (13) 

At large recoil region, the B c — > X c5 transition is dominated by the single gluon exchange 
as depicted FIGJ2 The expressions for F* c ^ Vc and A? f J/lp are listed in APPENDIX O 



F. The decay amplitudes 

The B c — > J/ipir, T] c ir decays are tree dominated within the framework of Operator Prod- 
uct Expansion, and without pollution from penguins and annihilation diagrams. In the 
perturbative QCD approach, the Feynman diagrams are shown in FIG]3l where (a) and (b) 
are factorizable topology, (c) and (d) are nonfactorizable topology. After a straightforward 
calculation using the modified perturbative QCD formalism Eq.(j3]), we obtain the decay 
amplitudes 

A{B-^X c ^-) = C ^V cb V: d £ A F1Gm , (14) 

v i=a,b,c,d 

where the CKM matrix elements V cb V* d = A\ 2 (l-A 2 /2-A 4 /8) + £>(A 8 ) with the Wolfenstein 
parameterization. The detailed expressions of ^FiGEb) are snown m APPENDIX [Dl From 
the expressions, we can clearly see that only the twist-2 distribution amplitude of the n 
meson contribute to the decay amplitudes. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The branching ratios in the B c meson rest frame can be written as: 

BK(B c ^X c - c it) = ^^-\A(B C ^X C - C n)\ 2 , (15) 



57T m 



B,- 



where the common momentum \p\ = (m 2 B — m| c J/2ms c , thejifetime and mass of the B c 
meson are rriB c = 6.276 ± 0.004 GeV and tb c = 0.46±0.07 ps [20(, respectively. Other input 
parameters are 

m c = 1.5 GeV, mj^ = 3096.916±0.011 MeV [20], fj^ = 405±14 MeV [18], 

m b = 4.20±g;5? GeV [20], m Vc = 2980.3±1.2 MeV [20], f Vc = 420±50 MeV [18], 

A = 0.814™ [20], A = 0.2257±ffi [20], f Be = 489±4 MeV [2JJ. 
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If not specified explicitly, we shall take their central values as the default input. 

Our numerical results a on the form factors F^ c ~* Vc and Aq\~^ J ^ are listed in TABLEjTJ 
where u and v are the parameters in the wave functions of Eqs. (lB7p - (1B10p and Eqs. (lB17p - 
flB20p . respectively. From the numbers in TABLElU we can see that 

(i) The form factors F and A),i,2 decrease with the increasing parameters u and v. The 
form factors of Fq and Ao,i are more sensitive to the parameters to or/and v than 
the form factor of A 2 . Uncertainties of form factors F and A 01 subjected to the 
charmonium wave function for Coulomb potential are larger than those for harmonic 
oscillator potential. Uncertainties of Fq and Ao t i related to the parameters us in our 
given range are about 16% ~ 20%, while those related to the parameters v in our given 
range are about 20% ~ 30%. In addition, the uncertainties of the decay constants Jb c , 
fj/rp and f Vc will bring ~ 0.8%, 3% and 12% uncertainties to the form factors F and 
A012, respectively. 



ii) The form factors 



27 



28 



29 



30 



lave been widely studied in the previous works [22|, |23j, |2J, |25|, [26, 



31 



22 



32]. There are large difference among the predictions in respect 



of various approaches. Compared with the previous results where F ( 







32J, our numerical results show that F ( 



B c — >r/ c 





A 



M l26|, l27j, |28|, [29|, [31 
^ > A^ J/tp . With some appropriate parameters, our results b on the form 



factors Fn and Aq^ are in agreement with those in the previous works 



22 



25 




26j, 



re form factors A 2 are smaller than those in the 



previous works |25l . 
there might be F c ^ Vc 



(or A»£ J/ *] 



321 ] . According to the "spectator quark" ansatz, 
or Aq^ 2 D *) ~ 0-6 by intuition. So 



IT'S— >_D 
^0 



maybe the results based on the three-point QCD sum rules 



23 



24| are small. 



Our numerical results on the amplitudes and branching ratios for B c — > J / i[)tt, r] c ir decays 



Here, we think that theoretical prediction on input parameters, such as lu and v, relies on our educated 
guesswork. All values within allowed ranges should be treated on an equal footing, irrespective of how 
close they are from the edges of the allowed range. For example, we cannot say that the probability of 
uj = 0.5 GeV is less than that of lu = 0.6 GeV, while the error means the usual one standard deviation 
in the form of j4j^ (such as the expression of ms c = 6.276 ± 0.004 GeV). So our numerical results had 
better to be given by a range to show the theoretical uncertainties, rather than the form of Aj^. 
For example, F^ Vc = 0.430 (0.464), A^ JH ' = 0.446 (0.470), A^ J/i) = 0.392 (0.427) for v = 0.80 
{lo = 1.60 GeV). 
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are listed in TABLEjTT] and IHII From the numbers in TABLEjTT] and 1111} we can see that 

(i) The contributions of the nonfactorizable topology ]FIG|3](c) and (d)] can provide large 
strong phases. The strong phases of the FIGJ3](c) topology 5 > 110°, while the strong 
phases of the FIGJ3] (d) topology 5 ^ —50°. The interferences between FIGJ3] (c) and 
(d) are destructive. The strong phases from nonfactorizable topology decrease with 
the increasing parameters u and v. They are free from the uncertainties of the decay 
constants /^ c , /jm and f Vc . The strong phases subjected to the charmonium wave 
function for Coulomb potential are larger than those for harmonic oscillator potential 
in our given ranges. 

(ii) The dominated contributions to the branching ratios come from the factorizable topol- 
ogy [FIG 13 (a) and (b)]. The ratio of amplitudes |^ F iG[3tc+d)/^FiG[3ta+b)l ~ f° r &c 
— > J/ipir decay, and about 2% ~ 3% for B c — > r] c n decay. The dominating amplitudes 
-4.Fic[2iab) and the branching ratios decrease with the increasing parameters u and v. 
Besides the large uncertainties from the parameter u and v, the uncertainties of the 
decay constants fj/^ and f Vt . will bring ~ 7% and ~ 24% uncertainties to the branching 
ratio for B c — > J/tpn and i] c tt decays, respectively. Considered the uncertainties from 
the input parameters, our results on the branching ratios are basically consistent with 



those in previous works 



31 



34 



Compared with the results in 



27 



|4l|, I42J, |43j (see the numbers in TABLEjIV]). 



2911 where small form factors are used (see the num- 



bers in TABLEjT]), we find that our predictions are large. If with the same factor 



factors, our results generally agree with those in [27|, |29j. The large predictions in [40] 
are obtained by the relations among the amplitudes under the quark diagram scheme, 
i.e. A(B c ^r) c ir) = V cb /V ub A(B^D s n) and A{B c ^J/^u) = V cb /V ub A{B^D s p). In- 
tuitively, the distribution amplitudes of the heavy quarkonia (such as B c , J/ip and 
rj c mesons) should be narrower than those of the "heavy-light" systems (such as B 
and D mesons). So, the superposition among the B c — J/tp(r] c ) — n systems might be 
less than those among the B — D s — ir(p) systems, i.e. there might be A{B c -^r] c 'K) < 
V cb /V ub A(B^-D s 7r) and A(B c ^J/^tt) < V cb /V ub A(B^D s pYlf this argument or/and 



assumption is true, then it is expected that the results in [40] would become smaller 
and be consistent with ours. 

(iii) The dominating amplitudes ^FicOab) an< ^ ^ ne branching ratios subjected to the char- 
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monium wave function for Coulomb potential are larger than those for harmonic oscil- 
lator potential. For a fixed value of parameter u or v, the relation of the branching ra- 
tios is B7Z(B c —>ri c iT) > BTZ(B c —> J/ijm). There are at least two reasons, one is that the 
phase spaces for B c — ► i] c 7r decay is larger than those for B c — > J j ipir decay, the other is 
that F^ v ° > A^ J '^ (see the numbers in TABLED. The relation of the branching 

Si- 



22 



27 



29 



37 



38 



40, 



41 



42 



ratios is in agreement with the previous predictions 
The signal of B c — ► J/ijj7r decay has been identified by the detectors at hadron collider 
Tevatron. It is eagerly expected that the signal of B c — > i] c 7r decay is at the near corner 
for Tevatron and LHC. 



IV. SUMMARY AND CONCLUSION 

In this paper, the B c — ► J/ip7T, r] c 7T decays are studied with the perturbative QCD ap- 
proach. It is found that the form factors Aq°{^ J and Fq c ^ Vc for the B c — > J/V>, r? c tran- 
sitions and the branching ratios for B c — > J/ipn and r] c 7c decays, they decrease with the 
increasing parameters uo and v, where uo and v are the parameters of the charmonium wave 
functions for Coulomb potential and harmonic oscillator potential, respectively. Therefore, 
the B c — > J/ipir, r] c TT decay modes provide good places to test quark potential models. In 
addition, the large uncertainties come from the uncertainties of the decay constants /ju 
and f Vc , which could be reduced greatly with the more accurate experimental measurements 
or/and better theoretical calculations. There are some other uncertainties not considered 
here, such as the power suppressed terms, the high order corrections, the effects of the final 
states interaction, the relativistic corrections to the wave functions, the model dependencies 
of the wave functions, and so on. They might be important in some cases (for example, 
the chirally-enhanced power corrections to the B — > nK decays are not much suppressed 
numerically.) and deserve the dedicated researches. So our results might be regarded as 
the estimations under the pQCD framework. One should not be too serious about these 
numbers. Anyway, the large branching ratios and the clear signals of the final states make 
the measurement of the interesting B c — > J/ipir, r\ c ii decays easily at the hadron colliders. 
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APPENDIX A: WAVE FUNCTIONS OF THE vr MESON 



The distribution amplitude for the twist-2 wave function and the distribution ampli 



tude (tP and cbt for the twist-3 wave functions are 



3 



_^6xx{l + 0.44C 2 3/2 (x -x) + 0.25C 4 3/2 (x -x)\, (Al) 
^j={l + 0.43C 2 1/2 (x -x) + 0.09Cl /2 (x - x)}, (A2) 



{c} /2 (z - x) + 0.55C 3 1/2 (x - x)}, 



2^/2N c 

with the decay constant f n = 130 MeV. The Gegenbauer polynomials are defined by 

1,_ , 



15 

T 1 



(A3) 



C{'\z) 


= z, 


c y\z) 


cy\z) 


= l(5z 3 -3z), 




cT{z) 


= 1(^-1), 


cT{z 



APPENDIX B: WAVE FUNCTIONS OF THE J/if> AND i] c MESONS 



The heavy quarkonium, such as cc, similar to diatomic molecules, might be amenable to 
a Born-Oppenheimer treatment c [l]. Following the prescription in 18, [33 1 , two forms of 
the wave functions corresponding to two different nonrelativistic potentials will be derived. 



1. wave functions for harmonic oscillator potential 

In the nuclear shell model, a more realistic description of the nucleons inside the atomic 
nucleus is given by the Woods-Saxon potential. The Schrodinger equation subjected to the 
Woods-Saxon potential cannot be solved analytically, and must be treated numerically, but 
the energy levels as well as other properties can be arrived at by approximating the model 
with a three-dimensional harmonic oscillator. The spectroscopy of the heavy quarkonium cc 
can be treated by this model. The quantum number n r L for the J/ip and r\ c mesons is IS, 

c The heavy quark-antiquark pair is bound by the gluon and light-quark clouds. The heavy quarks corre- 
spond to the nuclei in diatomic molecules. The gluon and light-quark fields correspond to the electrons, 
and provide adiabatic potentials [lj]. 
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where n r and L are the radial quantum number and the orbital angular momentum quantum 
number, respectively, (note : the energy spectrum of a three-dimensional harmonic oscillator 
is given by E nrL = {2{n r — 1) + L + §}o;.) The radial wave function of the corresponding 
Schrodinger state is given by 

4>n r L(r) = V>is(0 oc exp(-aV/2), 

where a 2 = m c oo/2, uj is the frequency of oscillations or the quantum of energy. 



:bi) 



Applying the Fourier transform, the state Eq.f lBlj) is replaced by the mapping represen- 
tation on the momentum space, 

■-P 
2a 2 



Employing the substitution ansatz |l8l. l33|: 
k i i , k 



d 3 r e ^ k i)i S (r) oc exp(— T ). 



(B2) 



( _ -\ m ° 

I OC OC } « 

v i 2 > 



ml + k 2 . 



xx 



(B3) 



where x — 1 — x, and x is the longitudinal momentum fraction of the valence quark of the 
meson, the wave function can be taken as 

(x — x) 2 m 2 + k\ 



i>is(k) -> ipis(x,k±) oc exp 



3a 2 xx 



(B4) 



Applying the Fourier transform to replace the transverse momentum k± with its conjugate 
variable 6, the 1 ^-oscillator wave function can be taken as 



ipis{x,b) ~ J d 2 k±e t ^ ± tpis(x,k±) oc xxexpj — 
The modified wave functions can be written as 

i>x ce (iS)(x,b) oc $ asy (x)exp{ - — xx[( 



-XX 
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V 2xx ) +UJ J J" 



X — X\ 2 



2xx 



) +u,V]}, 



(B5) 



(B6) 



with $ asy (x) being set to the asymptotic models of the corresponding twists for light mesons, 



which have been given in 
T] c mesons in Eq.((6]) and Eq.Q 

^M) " /W 
<t>%{x,b) 



331 ] . Therefore, we can obtain the wave functions of the J/ip and 



2y/2N c "'' 

fj/-4> 
2V2N C 



-xx 



A^xxexpj 

N^(x — x) expj xx ^ 



2xx 

m c _[ /x — x\2 

CO 



fvc AT v - f m c -17 

iv^xxexpj xx ^ 



CM) 



2^2N t 

fvc 



2xx 

X — X\ 2 



10' 



2^2N t 



ATS f m c - 

:iv„ exp<^ xx 



2xx 

X — X\ 2 

2xx 



(B7) 
(B8) 
(B9) 
(BIO) 
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where N c is the color number, N^, N^f are the normalization constants. All wave function 
in Eqs.( ]B7j) -( lB10l) are symmetric under x x and normalized : 

'dx <j)bHx,ti) = (Bll) 

'da; <p v > s (x,0) = (B12) 

The parameter u; « m^s) - mj/^ 1S ) « m^s) - ^ c (is) ~ 0.6 GeV. 



2. wave functions for Coulomb potential 

In the static QCD potential, the interactions between heavy quarkonium can be param- 
eterized and well described by a funnel shape Coulomb plus linear potential. At short dis- 
tances one-gluon-exchange leads to the Coulomb-like potential with a strength proportional 



to the QCD coupling constant a s |l| 

V(r) = -C F ^2 (B13) 
r 

where Cf = 4/3 is the SU(3) colour factor. 

The radial wave function of the corresponding Schrodinger state is given by (note : the 
principle quantum number n associated with Coulomb potential is given by n = (n r — 1) + 
L + l) 

4>n r L(r) = ^ts{r) oc exp(-g B r), (B14) 

where qs = Cpfi c a s is the Bohr momentum, /i c = m c /2 is the reduced mass of the c-quark. 
Analogous to the treatment for the case of harmonic oscillator discussed above, we can get 

^® K m 1 2N2 » (B15) 

(k 2 + q%) 2 

ip ls (x,b) oc v L c K 1 {m c byjl-4xx{l-v 2 )), (B16) 

1 — ixx(l — V 2 ) 



where the typical velocity of the quarks in charmonium v = qsj^c = 2a s /3 ~ 0.3 34|. The 
wave functions of the J/ ip and 7] c mesons can be written as 



l i .. M _ fm \ 1 (■'■■rfvij) 



b ) = V1=K I == / iKV 1 - 4 ^( X - « 2 )). ( B17 ) 
z v ZJV c 4/1 — Axx 1 — v 2 ) 



= vjl=K ) X ^^ K^m^l-Axxil-v 2 )), (B18) 
Zy/Zl\ c J\ _ ±xx(\ - V 2 ) 
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C (x, 6) = -fy=N? mcb KAmJoJl - 4xx(l - v 2 )), (B19) 

^ ; 2^ * ./l- 4*5(1 u cV 1 K J 



W = ^1 - 4**(1 - (B20) 



The normalization conditions are the same as those of Eq.( 1Bll[ ) and Eq.( 1B12l) . 



APPENDIX C: FORM FACTORS IN THE PERTURB ATIVE QCD APPROACH 



pBc-^Vc _ g nm 2 (j F I dxidx 2 / b 2 db 2 <j)B c {xi) 
Jo Jo 

x {E a (t a )a s {t a )H a (a,(3 a ,b 2 )[r Vc {2(l - x 2 ) - r^^fa, b 2 ) 

-((l-x 2 )-2r 6 )^ e (x 2 ,6 2 )' 
-E b (t b )a s (t b )H b (a, /3 b , b 2 ) [(rj^l - £i) + r c )#J c (x 2 , 6 2 ) 

-2r„ c ((l -xi) +r c )0^(x 2 ,6 2 )]} 



mft + mj/^ ^ = Cp ^ dxida , 2 f\ db2(j)Bc{xi) 

Jo Jo 



2m 



x {£ a (t a )a s (t a )# a (a, /3 a , 62) [(2 - x 2 - Ar b - x 2 r 2 ^(f)^(x 2 , b 2 ) 



+ (r&r^ - 2r^ + 4x 2 r v , + — )<%(x 2 , 62) 

-E b (t b )a s (t b )H b (a } /3 b , b 2 )(l + 2r c - 2xi + rj)0j;(x 2 , 6 2 )}, 



-4irm 2 Bc C F dxidx 2 b 2 db 2 4>b c (xi) 



™b c - mj/^ .B c ->j/ip 
-2mj/^ 2 ""Jo Jo 

x {E a (t a )a s (t a )H a (a,f3 aj b 2 ) ( - x 2 + x 2 rfy(j)fj,(x 2 , b 2 ) 

+ (?V> - 2fy H )^(^2, &2 

+E b (t b )a s (t b )H b (a,p b , b 2 )(l - 2 Xl )(l - rD^{x 2 ,b 2 )} 



where t 



0(6) 



max 



a|,V^|,l/6 2 ), ^a(b)(t) = e- s ^*), 
a = -m| c (a;i - x 2 )(a;i - r^x 2 ), 
Pa = -m|j(l - x 2 )(l - rjx 2 ) -r 2 b ] 
Pb = -m 2 Bc [{l-xx){rl-xx) -r 2 c ], 

S^(t) = s(x 2 p 2 , b 2 ) + s(x 2 p 2 , b 2 ) + 2 



* dfi 



1/62 /x 



■7 9 - 



(CI) 



(C2) 



(C3) 

(C4) 
(C5) 
(C6) 
(C7) 
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35]. The hard functions H are 



The quark anomalous dimension j q = — aJir. The explicit expression of s(Q,b) appearing 
in Sudakov form factor can be found in 

H a {a,(3,b) = 

H b (a,P,b) = 



K (b^) - K Q (by/JT) 

(3 — a 
K (b^) 







(C8) 
(C9) 



APPENDIX D: THE DECAY AMPLITUDES 



1. The amplitudes for B c — > J/ipir decay with the perturbative QCD approach 



/l POO 
^dx 1 dx 2 J ^b 2 db 2 (j) Bc (xi)E a (t a ) 

x a s (t a )C a (t a )H a (a, p^b^^r^l - x 2 ) - r 6 ]0j,(ar 2 , b 2 ) 

-[(l-x 2 )-2r b ]^(x 2 ,b 2 )}, (Dl) 

/l poo 
^dxidx 2 J Q b 2 db 2 cj) Bc (xi)E b (t b ) 

x a s (t b )C b (t b )H b (a,(3 b ,b 2 ){rl(l - Xl ) + r c }^(x 2 , 6 2 ), (D2) 

32tiCf f 1 f°° 

^figOc) = ^2W m Bc( l ~ r l) J o dxidx 2 da;3 J q b 2 db 2 b 3 db 3 Bc (xi)^(x 3 ) 

x E c (t c )a s (t c )C c {t c )H c (a, (3 C , b 2 , b 3 ){r^(xi - x 2 )cf)\{x 2 , b 2 ) 

+ (l-r})(xi-x 8 )^(x 2 ,fe)} l (D3) 

Aiclld) = ^=m Bc (l-rl)j dx^x^x^j ^b 2 db 2 b 3 db 3 (f) Bc (xi)(l) a w (x 3 ) 

x E d (t d )a s (t d )C d (t d )H d (a, (3 d , b 2 , h)\r^{x x - x 2 )^(x 2 , 6 2 ) 

+ [2{x 2 - - (x 2 - x 3 )(l - rj)]0j(a; 2 , b 2 )}, (D4) 



where t c(d) = max(J\a\, J\P c (d)\, 1/&2, 1/&3), -#c(d)(*) = e" 



-S^(t)-S^(t) 



Pc = -m 2 Bc (x 2 - xi)[{x 2 - xi) - (x 2 - x 3 )(l - rj)], (D5) 

/3 d = -m Bc (x 2 - xi)[(x 2 - xi) - (x 2 - x 3 )(l - rfy], (D6) 

/' du 

— 7g> ( D7 ) 
1/63 /i 

Ca(b) = Ci + C 2 /N C} C c (d) = C 2 , (D8) 
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H c{d) {a,P,b 2 ,h) = {e(b 2 -b 3 )K (b 2V ^)I (b 3V ^) + (b 2 ^b 3 )} 

x {e( +( 3)K (b 3 ^) + ^6(-(3)H^(b 3 ^)}, (D9) 

where C± j2 are the Wilson coefficients. The definitions of other parameters are the same as 
those in APPENDIX O 

2. The amplitudes for B c — > r] c ir decay with the perturbative QCD approach 



/l /*oo 
^dxida;2 J b 2 db 2 <f>B c (xi)E a (t a ) 

x a s {t a )C a {t a )H a {a,(5 a ,b 2 )[r Vc [2{\ - x 2 ) - r b ]<f>^(x 2 , b 2 ) 

-[(l-x 2 )-2r b }4>l(x 2 ,b 2 )}, (D10) 

/l yoo 
^dx 1 dx 2 J b 2 db 2 (j)B c (xi)E b (t b ) 

x a s {t b )C b {t b )H b (a, (3 b , b 2 ){[r* c (l - x x ) + r c ]<^ c (x 2) b 2 ) 

-2r Vc [(l-x l ) + r c }<p s Vc (x 2 ,b 2 )}, (Dll) 

— i327rCp r 1 r 00 

^figHc) = y^i\r m BcO-- r D J dxidx 2 dx 3 J o b 2 db 2 b 3 db 3 (j)B c ( x i)K( x s) 

x E c (t c )a s (t c )C c (t c )H c (a, (3 C , b 2l b 3 )\r Vc {xx ~ x 2 )^(x 2 , h) 

- [(1 - r-y (x! - x 3 ) + 2rJ> 1 - x 2 )]0;(x 2 , 6a)}, (D12) 
^FicUd) = ^2N~ ~ r i) J dx 1 dx 2 dx 3 J J) 2 db 2 b 3 db 3 (j) Bc {xi)(j)%{x 3 ) 

x E d (t d )a s (t (i )C ( i(t d )i^ d (a,/5 d ,62,&3){^ c (xi - x 2 )0* o (x 2 , 62) 

+ [2(x 2 - xx) - (x 2 - x 3 )(l - <)]</>; (x 2 , 6a)}. (D13) 
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TABLE I: Form factors of Rf c ^ c and A. 





D frt 

T?£>c >r lc 


^0 


. R —> T/ih 
c J i ™ 




oj = 0.5 GeV 


0.790 


0.775 


0.671 


0.469 


uj = 0.6 GeV 


0.741 


0.730 


0.636 


0.454 


oj = 0.7 GeV 


0.698 


0.690 


0.605 


0.440 


uj = 0.8 GeV 


0.660 


0.655 


0.578 


0.427 


v = 0.25 


0.903 


0.891 


0.712 


0.363 


v = 0.30 


0.824 


0.819 


0.664 


0.364 


v = 0.35 


0.760 


0.759 


0.624 


0.361 


w = 0.40 


0.705 


0.708 


0.589 


0.356 


[22] 


0.170 ~ 0.687 


0.156 ~ 0.684 


0.156 ~ 0.745 


0.156 ~ 0.862 


[23] 


0.20±0.02 


0.26±0.07 


0.27±0.03 


0.28±0.09 


[24J 


0.23±0.01 


0.21±0.03 


0.21±0.02 


0.22±0.02 


[25] 


0.66 


0.60 


0.63 


0.69 


[26] 


0.76 


0.69 


0.68 


0.66 


[2D 


0.47 


0.40 


0.50 


0.73 


[28J 


U.Dl 


U.00 


U.00 


U.04 


[29] 


0.49 


0.45 


0.49 


0.56 


[30] 


0.87 


0.27 


0.75 


1.69 


[31] 




^7+ - 01 

u -°'-0.02 


^+ - 01 

u - oo -0.03 


u - oi -0.04 


[32] 


n si +0.01+0.03 

u - ui -0.02-0.04 


0.53±0.01 


^n +0 - 01 

u - ou -0.02 


n 44+ - 02 
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TABLE II: Amplitudes and branching ratio for B c — > J /ipir decay, where S is the strong phase. 





AlGUa) 




^figHc) [<*] 


A 


^IGUd) I s . 




BK{B c -^J/ipir) 






0.5 GeV 


1.359 


1.831 


-0.115 + i0.269 [113°] 


+0.132 


- i0.285 [ 


-65°] 


1.913x 10" 


_3 


U) 




0.6 GeV 


1.235 


1.767 


-0.103 + i0.267 [111°] 


+0.115 


- i0.285 [ 


-68°] 


1.689x 10" 


_g 


U 




0.7 GeV 


1.133 


1.704 


-0.093 + zO.261 [110°] 


+0.101 


- i0.279 [ 


-70°] 


1.506x10" 


-3 


uJ 




0.8 GeV 


1.049 


1.643 


-0.086 + i0.253 [109°] 


+0.090 


- &0.271 [ 


-72°] 


1.352x10" 


-3 




V 


= 0.25 


1.941 


1.738 


-0.140 + f0.150 [133°] 


+0.157 


- i0.192 [ 


-51°] 


2.542x10" 


-3 




V 


= 0.30 


1.692 


1.686 


-0.130 + i0.162 [129°] 


+0.144 


- i0.203 [ 


-55°] 


2.140x10" 


-3 




V 


= 0.35 


1.498 


1.632 


-0.121 + i0.169 [126°] 


+0.130 


- i0.209 [ 


-58°] 


1.834x10" 


-3 




V 


= 0.40 


1.342 


1.576 


-0.111 + i0.174 [123°] 


+0.118 


- t0.212 [ 


-61°] 


1.591x10" 


-3 



TABLE III: Amplitudes and branching ratio for B c — > t] c tt decay, where 5 is the strong phase. 





-^FIGHa) 


-^FicHb) 


-^FIGHc) 


[S] 


-^FIGHd) 


[S] 


BTZ{B c ^ Vc 


7T 






0.5 GeV 


1.490 


1.828 


-0.112 + i0. 211 


[118°] 


+0.127 


- i0.300 


[-67°] 


2.117x10" 


-3 


uJ 




0.6 GeV 


1.353 


1.756 


-0.098 + i0.214 


[115°] 


+0.112 


- i0.299 


[-69°] 


1.858x10" 


-3 


LU 




0.7 GeV 


1.242 


1.685 


-0.087 + i0.212 


[112°] 


+0.100 


- i0.292 


[-71°] 


1.646x10" 


-3 


UJ 




0.8 GeV 


1.149 


1.617 


-0.078 + i0.207 


[111°] 


+0.089 


- i0.283 


[-73°] 


1.470x10" 


-3 




V 


= 0.25 


2.140 


1.665 


-0.131 + i0. 104 


[141°] 


+0.154 


-i0.189 


[-51°] 


2.792x10" 


-3 




V 


= 0.30 


1.864 


1.608 


-0.121 + i0. 119 


[135°] 


+0.141 


- i0.201 


[-55°] 


2.323x10" 


-3 




V 


= 0.35 


1.649 


1.548 


-0.112 + i0. 129 


[131°] 


+0.129 


- f0.208 


[-58°] 


1.970x10" 


-3 




V 


= 0.40 


1.476 


1.488 


-0.102 + i0.137 


[127°] 


+0.118 


-i0.212 


[-61°] 


1.692x10" 


-3 



TABLE IV: Branching ratio for B c — > J/ifiir, r] c ir decay in previous works (in the unit of 10 3 ). 



BTl{B c - 




0.13~1.55 [22] 
0.26 [39] 


0.85 [27] 
9.30 [40] 


0.94 [29] 
2.00 [41] 


1.1 1-2. 16 [36] 
1.16~1.34 [42] 


2.30 [37] 1.8 [38] 
1.90 [43] 


BU{B C - 


>J/V>tt) 


0.02~0.34 [22] 
1.30 [39] 


0.61 [27] 
4.50 [40] 


0.76 [29] 
1.30 [41] 


2-o±8;?±8;? [21] 

1.08~1.24 [42] 


2.19 [37] 1.7 [38] 
1.70 [43] 
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"P3 = Pi ~ P2 




(a) (b) 
FIG. 2: Feynman diagrams contributing to the B c — ► X cg form factors. 




(a) (b) (c) (d) 

FIG. 3: Feynman diagrams for B~ — ► X cg 7r~ decays. 
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